Variational approach for the two-dimensional trapped Bose Einstein condensate 
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We construct a many-body Gaussian variational approach for the two-dimensional trapped Bose 
gas in the condensate phase. Interaction between particles is modelized by a generalized pseudo- 
potential of zero range that allows recovering perturbative results in the ultra-dilute limit, while 
back action of non-condensate particles on the condensate part is taken into account for higher 
density. As an application, we derive the equation of state and solve stability problems encountered 
in similar mean-field formalisms for a single vortex configuration. 
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Both for fundamental reasons and for eventual fu- 
ture applications, ultra cold atomic gases trapped in re- 
duced geometries has become a challenging field. Two- 
dimensional (2D) configurations are especially interesting 
for their various expected features including the interplay 
between the Bose Einstein and the Kostcrlitz-Thouless 
transitions, the concept of tunable BEC [l| and in the ro- 
tating case the possible realization of a gas with fractional 
statistics quasi-2D BEC has been produced using 

polarized hydrogen atoms adsorbed on a helium film [3j, 
with clear evidence of the quasi-condensate phase. Sev- 
eral experimental groups 0, IE IE are now interested 
in quasi-2D configurations using optical and magneto- 
optical traps which are expected to be more flexible for 
varying physical parameters such as the atomic density 
or the interaction strength. On the theoretical side, the 
general formalisms deeply differ from usual many body 
treatments in that they take into account both thermal 
phase fluctuations and the specific scattering properties 

liinananina. 

In this Letter we use the Gaussian variational approach 
of the many-body problem to study the condensate phase 
for the 2D BEC We show that describing interac- 

tions between atoms with the so-called A-potential 0] 
is a very simple way to construct a gapless and consis- 
tent theory. In 3D configurations, this formalism has 
been already developed [15| and allows to recover in a 
purely variational way the gapless HFB-G2 prescription 
[l6j . In the 2D regime, this way of taking into account 
the back action of non-condensed particles on the conden- 
sate field has not been yet implemented and brings new 
features. The Letter is composed of three parts. In the 
first part, we describe the modehzation of the system and 
recall some basics of the two-body problem in 2D using 
the zero range approximation. In the second part, we 
develop the Gaussian variational approach (called also 
Hartree Fock Bogolubov or HFB) for the 2D trapped 
Bose gas and show the links with the HFB-G2 scheme. 
In the last part we apply the formalism in two distinct sit- 
uations. We first calculate the Equation Of State (EOS) 
at zero temperature in the homogeneous limit. This first 
result generalizes the known perturbative EOS. Finally, 
as a test of the formalism we study the thermodynami- 



cal instability of a single vortex in a non rotating trap. 
We conclude that unlike unexpected results obtained in 
similar but non variational approaches |17| . the system 
remains unstable at finite temperatures. 

We consider N bosons of mass m in a very anisotropic 
trap characterized by a high axial frequency tu z and small 
longitudinal frequencies oj x = u> y = lul <C w z . In the 
following, we suppose that the system is in a quasi two- 
dimensional geometry. This situation is achieved when 
the typical energy (e) of a two body scattering process 
is sufficiently low as compared to the axial frequency 
|e| <C hto z . In this regime, the system is somehow frozen 
along the tight direction in the ground state of the axial 
trapping potential and we can restrict our study to the 
longitudinal xy-plane. Before going into the approximate 
treatment of the many body problem, let us recall some 
basics of the two-body case in a planar trap (u>l = 0). 
We note the spatial coordinates in the plane (fi, fb), the 
relative and center of mass coordinates are respectively 

r = ?\ — r>. 2 and R = 1 — 2 . In the zero range ap- 
proach, the two-body wave function ip2 obeys the free 
Schrodinger equation for r =/= 0, while the interaction is 
taken into account by imposing the contact condition : 

Mn,r 2 ) = A(R)]n ( — ) for r^O , (1) 
\a.2D J 

where a,2D (the "2D scattering length" ) is obtained from 
the usual scattering length a^o and from the character- 
istic length of the trap a z using the relation 

a 2D = 2.092 a z exp (- ^| . (2) 

Note that the short range part of the wave function is 
not described by relation §| which gives only the 
correct asymptotic behavior for r 3> a z . However, the 
zero range approach is justified when the mean inter- 
particle distance is much larger than the axial width a z — 

In term of the 2D atomic density n, this gives 

na 2 z <C 1. Another and equivalent way to implement the 
zero-range approach is to solve the Schrodinger equation 
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using the A-potential [15| defined by: 



(i?,r|l/ A |^ 2 )=.g A( 5( 2 )(r)^ A 



R+t;) 



(3) 



In this potential, g\ is the two-body T-matrix evaluated 
at energy ~h 2 A 2 /m : 



1 



m ln(qAa2i)) 



(4) 



where q = e 7 /2, 7 is the Euler's constant and TZa is a 
regularizing operator: 



lim {1 

!-->0 



ln(qAr)r9 r } 



(5) 



By construction, the observables are independent of the 
parameter A which can be considered as a free field of 
the exact theory, function of the center of mass coordi- 
nates R 0| . However, we recall that some approximate 
treatments are A-dependent. In this type of formalism, 
A appears as a new degree of freedom which is a way of 
improving the results. Taking for example the first order 
Born approximation in the two-body case, one can easily 
check that the exact scattering amplitude at energy e is 
obtained using an uniform A = —i^/me/Ti. In this case, 
A is directly linked to the energy of the scattering process 
and has a straightforward physical meaning. 

Now, we implement the variational Gaussian approach 
in the case of a 2D condensate in the external poten- 
tial V^ x t = mui using the A-potential (|3J. In 
Ref. [13 , it has been already shown that this type of ap- 
proximation is A-dependent. We explain in the follow- 
ing the criterion we use for a systematic determination 
of the field A(R). In order to simplify the discussion, 
we break the U(l) symmetry by splitting the atomic 
field ip into a classical field $ and a quantum fluctuation 
cj) — tp — $ The HFB scheme consists in choosing 
as a trial density operator the exponential of the most 
general quadratic form of the fields {(f), (ft}. Diagonaliza- 
tion of the quadratic form is obtained using a Bogolubov 
transformation : 



W = E {bnUn(R) + Hy n (R)} 



(6) 



where {b n , b^} are the usual creation and annihilation 
bosonic operators for a quasi-particle of quantum num- 
ber n, while the amplitudes {u n ,v n } verify the modal 
equations : 



-—A + V; xt + TiSii - 
2m 



— A + V cxt + hS* n -n 
2m 



u n + hT, 12 v n = hu n u n 

(7) 

v n + Ti£* 12 u n = -huj n v n . 



From the variational principle 0] we deduce the self- 
energies : 



nEn(R) = 2g A (^ (RWM 



hY, 12 (R) = gA^A 



(8) 
(9) 



and also a generalized Gross-Pitaevskii equation for the 
classical field 



2m 



A + Vcxt + 2gA?i - /i 



$ + ?i£i2$* =0 



(10) 



where n is the non-condensed density. Let us insist on 
the fact that the modal equations J7J) do not coincide 
with the ones appearing in the linearized time-dependent 
treatment (the so-called RPA formalism; see for example 
Ref. 01)- a consequence, while the condensate mode 
associated to a phase change in the condensate wave func- 
tion is a solution of the time dependent equations with 
a zero eigen-frequency, it is not in general solution of 
Eqs.JJJ). In the homogeneous case, this leads to the pres- 
ence of a gap in the spectrum defined through Eqs.Q), 
and in the inhomogeneous case, to the presence of an 
unphysical mode of non-zero eigen-frequency. Unfortu- 
nately, the existence of a spurious energy scale introduced 
by this mode prevents from a consistent description of 
the thermodynamical properties of the system. We show 
in the following that this problem can be solved with- 
out modifying the equations deduced from the variational 
principle. For that purpose, we recall that Eq.JHJ corre- 
sponds to a first order Born approximation so that unlike 
an exact approach, HFB is explicitly A-dependent. The 
key point is then to choose a specific realization of the 
field A, noted in the following A*, such that the conden- 
sate mode: 



u (R) = $(R) ; vo(R) = -$*(R) 



(11) 



solution of the time dependent equations, is also solution 
of Eqs.Q with a zero frequenc y. T his requirement leads 
to the condition = 5a* < E >2 |2j, so that at each point 
R in the condensate, A*(R) is an implicit solution of the 
equation |23|: 



Ha* 



${R- r ~)^R+ r ~)) 







(12) 



In the limit 03.0 *C a z where interatomic collisions are 
weakly modified by the transverse confinement, Ea. H12|) 
leads at finite temperature to the HFB-G2 prescription 
[lr| introduced in the 3D configuration : 



A-kU a 3D 



ma 2 



$2 



(13) 



where kt(R) = Z^2 n ^9nU n (R)v^(R) is the thermal con- 
tribution of the pairing field and g n is the usual Bose 
occupation factor for the mode of energy hu) n . 
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FIG. 1: Equation Of State of the 2D homogeneous Bose 
gas. Dotted line: EOS deduced from Schick's formula /i = 
47rfo 2 n/ (m lnfwaor->'0|24|. Full line: full variation al ap proach 
(HFB). Dashed dotted line: Perturbative EOS fl FOf . 

As a first application, we consider the homogeneous 
medium at zero temperature. From Eas. (|7HUI12[l . we 
deduce a standard Bogolubov spectrum and a self- 
consistent equation for A*: 

^-=.9A*|$| 2 , (14) 
m 

giving directly the energy at which the T-matrix has to 
be evaluated for having a gapless spectrum. From fi = 
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FIG. 2: Typical evolution of the anomalous mode energy of 
a single vortex in a non rotating trap as a function of the 
temperature. The presence of a negative energy mode in the 
spectrum leads to the so-called thermodynamical instability of 
the single vortex |26jl . Dashed line: the coupling constant used 
in the gapless equations is arbitrarily evaluated using the two- 
body T-matrix at energy — 2/ii oc . Full line: full variational 
scheme; at finite temperature, the anomalous mode energy is 
not stabilized. 



(/A* C$1 2 + 2n) and (jT4*|) we deduce the EOS shown in 
Fig-©- In the ultra-dilute limit na\ D < 1, we recover 
the exact perturbative result 0,0: 

n= T^ ln {e^ m ,alJ ■ (15) 

For increasing values of the 2D atomic density, the pre- 
dicted EOS deviates from Ea. (|15(l (for a value of the 2D 
gas parameter greater than 3 10 -5 , the relative difference 
between the two EOS is more than 5%). Note that sim- 
ilarly to the 3D situation, there is an upper bound for 
the values of the density that can be studied 15], here 
at T — 0, n < n crit = exp(— 2j)/ira,2 D . At n = 7i. C ritj 
/i — > oo and $ — > 0, this gives a fundamental limit of 
applicability of the present formalism. In 2D traps, we 
have checked that at zero temperature, for static con- 
figurations, the full variational approach gives the same 
results as the LDA. In situations where the density varies 
rapidly, for example in presence of vortices, or also at fi- 
nite temperature not only Eas. (|7ll()fl but also Ea. (fT^I 
have to be solved numerically. We consider then as a last 
example, the case of a single vortex in a non rotating trap 
at finite temperature |25| . This configuration is especially 
interesting as a test of the present formalism. Indeed, a 
previous self-consistent but not variational approach |l7j ] 
has led to the conclusion that a vortex which is thermo- 
dynamically unstable at vanishing temperatures, could 
be stabilized at finite temperature. We recall that this 
instability is due to the existence of a core localized state 
having a negative energy (huj a ) when the trap rotation 
frequency is zero |2r| . In 3D rotating BEC experiments, 
this type of instability which causes the spiraling of the 
vortex core out of the condensate has been observed [2l} , 
and we expect a similar behavior in 2D. We then consider 
an off-equilibrium situation where the trap is at rest and 




FIG. 3: Typical density profile of a single vortex at finite tem- 
perature (fcsT = 20Huil), obtained using the full variational 
approach when the anomalous mode is not thermalized. In- 
set: mean-field coupling constant in harmonic oscillator units, 
obtained by solving Eg. 11211 . 
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all the modes except the anomalous core localized state 
are thermalized. In a first step, we have checked that ac- 
cordingly to Ref.^3], stabilization occurs when the cou- 
pling constant g\ is arbitrarily set to a constant value 
in the "gapless" HFB equations. In a second step, we 
have used a more realistic prescription by assuming that 
the T-matrix entering the "gapless" equations is taken 
at energy —h 2 A 2 /m = — 2/ii oc (/xjoc is the local chemi- 
cal potential). As shown in Fig[3 this non- variational 
approach predicts also an unphysical stabilization at suf- 
ficiently high temperature. On the contrary, using the 
full variational scheme where Ea. (|12[l is solved numeri- 
cally, the anomalous mode is never stabilized (see Fig|2J). 
In Fig.®, we show a typical density profile and the as- 
sociated mean-field coupling constant qa* obtained at 
finite temperature as a function of the distance from 
the core (we have used the notation — \J n / 'raw l) ■ 
As explained in Ref.^Tj; t ne unexpected stabilization in 
the non-variational approach is linked to the mean-field 
term g\h in Eg. (|lUfl which acts as a pinning potential. 
In the formalism presented here, the low energy modes 



which have a very different structure with respect to the 
homogeneous case, crucially determine the behavior of 
the coupling constant g\* through Ea. lfT^f . As a conse- 
quence, the homogeneous result Ea. l|14[) cannot be ap- 
plied: the lowest modes which become widely occupied 
with increasing temperatures, shape the structure of <?a* . 
For example in Fig.®, the coupling constant dives at the 
origin (recall that the structure at the vicinity of the core 
plays a central role in the pinning process) and presents 
a buttress at the edge of the condensate due to the pop- 
ulation of the surface modes. 

In conclusion, we have presented an approach for the 
2D trapped BEC which takes into account the back ac- 
tion of the non-condensed fraction on the condensate field 
at finite temperature. This formalism includes known 
results on the 2D BEC and solves stability problems en- 
countered in similar but non variational approaches. 
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